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PRELIMINARY  REPORT  OF  NUMERICAL  SIMULATIONS  OF  INTERMEDIATE  WAVE¬ 
LENGTH  ExB  GRADIENT  DRIFT  INSTABILITY  IN  IONOSPHERIC  PLASMA  CLOUDS 

INTRODUCTION 

Barium  cloud  injection  experiments  [Rosenberg ,  1971;  Davis  et  alM 
1974]  have  provided  much  data  concerning  the  dynamical  evolution  of 
plasma  clouds  in  the  F  region  ionosphere.  The  characteristic  steepen¬ 
ing,  elongation,  and  striation  of  barium  plasma  clouds  have  been  ex¬ 
plained  by  applying  the  ExB  gradient-drift  instability  [Simon,  1963] 
to  plasma  cloud  geometries  [Haerendel  et  al . ,  1967;  Linson  and  Workman, 
1970;  Volk  and  Haerendel,  1971;  Perkins  et  al . ,  1973].  Numerical 
solution  of  the  fundamental  fluid  equations  modelling  the  ExB  insta¬ 
bility  in  plasma  clouds  have  reproduced  not  only  the  large  scale  gross 
observational  features  of  plasma  cloud  evolution  [Zabusky  et  al.,  1973; 
Lloyd  and  Haerendel,  1973;  Goldman  et  al.,  1974;  Doles  et  al . ,  1976; 
Qssakow  et  al. ,  1975,  1977]  but  also  the  spatial  power  spectra  [Scanna- 
pieco  et  al . ,  1976]  and  minimum  scale  sizes  of  plasma  cloud  striations 
[McDonald  et  al.,  1978].  However,  the  goal  of  these  numerical  simula¬ 
tions  was  to  model  the  dynamics  of  the  entire  plasma  cloud  and  wave¬ 
lengths  less  than  about  a  kilometer  were  not  accurately  resolved.  But 
barium  cloud  experiments  sponsored  by  the  Defense  Nuclear  Agency  have 
shown  [Baker  and  Ulwick,  1978;  Kelley  et  al. ,  1979]  that  the  power 
spectra  of  the  plasma  cloud  striations  extends  to  wavelengths  on  the 
order  of  meters  to  tens  of  meters.  It  is  of  interest  to  numerically 
model  these  intermediate  wavelength  (100-1000  m)  irregularities  in  order 

to  compare  with  and  supplement  experimental  observations. 

Note:  Manuscript  submitted  October  5,  1979. 
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In  this  preliminary  report  we  present  a  two-dimensional  numerical 
simulation  of  the  intermediate  wavelength  Ex_B  gradient-drift  insta¬ 
bility  applicable  to  local  unstable  regions  of  large  F  region  barium 
plasma  clouds.  The  results  of  these  simulations  will  be  shown  to  be 
consistent  with  experimental  and  theoretical  studies  of  striated  barium 
clouds  and,  in  many  respects,  similar  to  recent  numerical  simulations 
fKeskinen  et  al..  1979]  of  the  intermediate  wavelength  collisional 
Fayleigh-Taylor  instability  in  local  regions  of  the  bottomside  of  the 
equatorial  F  region  ionosphere. 

MODEL  EQUATIONS 

We  wish  to  model  two-dimensional  Ex.B  gradient-drift  processes  in 
local  unstable  regions  of  large  ionospheric  F  region  plasma  clouds. 

The  restriction  to  large  clouds  (large  Pedersen  conductivity  compared 
with  that  of  the  background  ionosphere)  permits  both  the  neglect  of 
the  cloud  interaction  with  the  background  ionosphere  (second  level) 
and  variations  of  cloud  density  and  potential  along  the  magnetic 
field  lines.  For  wavelengths  greater  than  the  ion-gyroradius  (approxi¬ 
mately  10  m  for  Ba+  in  the  twilight  F  region)  a  fluid  description  can 
be  used  which  equations  have  been  given  previously  [Perkins  et  al., 
1973;  McDonald  et  al. ,  1978;  Chaturvedi  and  Ossakow,  1979]. 

By  adopting  a  Cartesian  coordinate  system  (x,y,z)  with  the  geo¬ 
magnetic  field  in  the  z-direction,  the  ambient  electric  field  E^ 
along  the  x-axis,  and  ignoring  to  lowest  order  electron  and  ion 
inertia,  we  can  write  after  transforming  to  the  cE  xB/B2  frame 
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9n 

3t 


(1) 


c  „  -  „  2cT 

-  -  Vcp^z-Vn  =  —  V2n 

e 


T 

V  •  nVcpi  +  —  V^n  =  •  Vn  (2) 

where  n(x,y,t)  is  the  plasma  cloud  ion  density,  cpi(x,y,t)  is  the 

electrostatic  potential  of  the  plasma  cloud,  K  ^  =  K  ^  +  K  .  \ 

e  en  ei 

K  =  ft  /v  ,  with  v  ,  v  .  the  electron  collision  frequencies  with 
ea  e  eci  en  ei 

neutrals  and  ions,  respectively,  and  =  T\  =  T  is  the  temperature. 

All  other  symbols  retain  their  conventional  meaning.  By  linearizing 

(1)  and  (2)  and  expressing  n  =  n^  +  nj  with  n^,  a  exp 

[i(k  x  +  k  y)  +  Vt  t],k*B  =  0,kL  >>1  we  find 
x  y  k.  —  — 


Y,  =  (cE  /BL)  (k  /k)2  -  (v  +  V  . )  (k2C  2)/n  n. 
k  o  x  en  ei  s  e  l 


where  L  1  =  (1/n  )  (3n  /3y) ,  k2  =  k  z  +  k  %  C  z  =  T/m. .  We  note  that 
o  oy  x  y  s  i 

the  growth  rate  y^  maximizes  for  modes  perpendicular  to  the  density 
gradient  (k  =  k^)  while  the  modes  parallel  to  the  density  gradient 
(k  =  ky)  are  damped  by  cross-field  diffusion. 

Numerical  Simulations 


By  defining  n*  (x,y)  =  n (x,y) /nQ (y) ,  cpj  (x,y)  =  cp1(x,y)/BL, 
x*  =  x/L,  y*  =  y/L,  t*  =  ct/L,  where  nQ(y)  will  be  defined  later, 
equations  (1)  and  (2)  can  be  written  in  dimensionless  form  as  follows 
(after  dropping  primes  for  clarity) 
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WLl  ,  LA  V  1.1  W  v 

3y  +  ^  =  '  ^2 


I  *IL 

n  3x  (4) 


where  =  2T/eBK^L,  f32  =  E  /B,  £3  =  T/eBL  are  dimensionless  constants. 

Equations  (3)  and  (4)  were  solved  numerically.  The  normalized 

plasma  cloud  density  n/n^  in  (3)  was  advanced  in  time  using  a  variable 

time  step  flux-corrected  predictor-corrector  scheme  [Zalesak,  1979] 

which  is  basically  second  order  in  time  and  fourth  order  in  space. 

Equation  (4)  was  solved  for  the  cloud  potential  cp^(x,y)  at  each  time 

step  using  a  regridded  Chebychev  relaxation  method  [McDonald ,  1977]. 

The  results  to  be  presented  here  were  obtained  using  mesh  of  64  x  64 

grid  points  with  constant  spacing  Ax  =  Ay  =  15  m.  Periodic  boundary 

conditions  were  imposed  on  n/n^  and  cp  ^  in  both  the  x  and  y  directions. 

The  parameters  used  in  these  simulations  are  typical  of  200  km  barium 

releases:  T  =  103  °K,  ^  e~  9  x  10^  sec  1 ,  (Ba"*")  =  36  sec  1 ,  B  = 

0.5  G,  E  =  2mV/m,  v  .  =  2  x  103  sec  v  =  150  sec  1  [Banks  and 
o  ei  en  - 

Kockarts ,  1973;  McDonald  et  al.,  1978]. 
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Three  different  computer  runs  were  made  using  different  values 

of  the  equilibrium  plasma  cloud  density  gradient  scale  length  L  =  3, 

6,  10  km  where  nQ(y)  =  Nq  (1  +  y/L) ,  constant,  is  a  steady  state 

( 3 / 3 1  =  0)  solution  of  (1)  and  (2).  The  simulations  were  initialized 

with  a  two-dimensional  perturbation  of  the  form  [Chaturvedi  and 

Ossakow,  1979]  6n(x,y,t  =  0)/n  =  A  sin  k  y  cosk  x  +  A2  osin  2k  y 

o  1  >  1  y  x  ,  y 

with  k  =  k  =  2tt/960  m  with  A^  1  =  2  x  10  5  and  A2  q  =  2  x  10  &  (where 
x  y  9  9 

6n  =  n-n^).  We  will  now  present  the  important  nonlinear  results  of 
these  simulations. 

Fig.  la  gives  an  isodensity  contour  plot  of  6n((x,y)/nQ  in  the 
x-y  plane  at  t  =  0  for  L  =  3  km.  The  initial  contours  describe  a 
sequence  of  enhancements  (6n/n^  >  0)  and  depletions  (6n/n^  <  0) 
arranged  in  a  checkerboard  fashion.  Fig.  lb  and  lc  show  the  evolution 
of  6n/n at  t  =  1500  and  1800  sec,  respectively.  The  density  fluctua¬ 
tion  contours  of  6n/n^  at  t  =  2000  sec  are  displayed  in  Fig.  Id  where 
some  elongation  in  the  y-direction  (ExB)  and  steepening  can  be  seen. 
Similar  contour  development  was  also  observed  for  the  other  two  plasma 
cloud  density  gradient  scale  lengths  L  =  6,  10  km  but  on  longer  time 
scales . 

Figs.  2a-b  show  representative  one-dimensional  power  spectra  both 
parallel  P(k  )  to  the  ExB  drift  and  perpendicular  P(k  )  in  the  non- 

y  x 

linear  late  time  regime  for  L  =  6  km.  These  power  spectra  are  defined 
as  follows 
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P(k  )  =  /  dk  I  6n(k  ,k  )/n  \  2 

x  /  y  x  y  o' 

P(k  )  =  /  dk  I  6n(k  ,k  )/n  \2 

y  /  x  1  x  y  o 1 

In  the  late  time  nonlinear  state  for  the  three  scale  lengths  L  =  3,  6, 
-n  -n 

10  km,  P(k  )  k  X  and  P(k  )  k  ^  with  n  =**  2-3  and  n  ^  2-2.5 
xx  y  y  x  y 

for  wavelengths  X  =  70-960  m.  These  spectral  indices  are  consistent 
with  recent  in  situ  experimental  measurements  [Baker  and  Ulwick,  1978; 
Kelley  et  al.,  1979]  and  analytical  studies  [Chaturvedi  and  Ossakow, 
1979]  of  barium  cloud  striations. 

Finally,  Fig.  3  presents  the  time  history  of  modes  Aj  \  and  A2  o 
for  L  =  10  km.  Initially  the  A2  o  mode  along  the  .ExB  (density  gradi- 
ent)  direction  damps  as  given  by  linear  theory.  As  the  linearly  un¬ 
stable  Aj  i  increases  it  nonlinearly  triggers  the  growth  of  A2  q •  At 
late  times  these  modes  arrange  themselves  in  approximate  agreement  with 
the  nonlinear  amplitudes  predicted  by  Chaturvedi  and  Ossakow  [1979]. 
Similar  time  development  of  A^  i  and  A2  0  are  observed  for  L  =  3,  6  km. 

j  > 

SUMMARY 

We  have  performed  preliminary  numerical  simulations  of  the  inter¬ 
mediate  wavelength  _ExB  gradient-drift  instability  in  local  unstable 
regions  of  large  F  region  ionospheric  plasma  clouds  using  parameters 
typical  of  200  km  barium  releases.  For  barium  cloud  density  gradient 
scale  lengths  L  =3,  6,  10  km  we  find  that:  (1)  linear  unstable 
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modes  saturate  by  nonlinear  excitation  of  linearly  damped  modes  along  the 
y  direction  parallel  to  ExB  drift  and  the  initial  cloud  density  gradient; 
(2)  in  the  nonlinear  well  developed  state  elongation  along  the  _ExB  di¬ 
rection  is  seen;  (3)  the  one-dimensional  power  spectra  P(k  )  «  k  nx  and 
-n  xx 

P(k  )  «  k  ^  with  n  ~  2-3  and  n  -  2-2.5  for  wavelengths  X  =  70-960  m. 
y  y  x  y 

These  results  are  consistent  with  recent  barium  cloud  experimental  measure¬ 
ments  [Baker  and  Ulwick,  1978;  Kelley  et  al. ,  1979];  support  the  nonlinear 
analytical  work  of  Chaturvedi  and  Ossakow  [1979];  and  show  furthermore 
the  similarity  between  the  ExB  gradient  drift  instability  in  ionospheric 
plasma  clouds  and  the  collisional  Rayleigh-Taylor  instability  in  the 
equatorial  F  region  ionosphere. 

Future  studies  are  planned  which  include  variation  of  the  para¬ 
meters  used  in  these  simulations,  addition  of  inertial  effects,  and 
inclusion  of  the  cloud  interaction  with  the  background  ionosphere 
(second  level). 
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Fig.  3  —  Time  history  of  mode  amplitudes  A2io  (along  initial  cloud  density  gradient) 
and  Ai  i  for  L  =  10  km  where  Am>n  =  Amkoy>nkox  and  kox  =  koy  =  27T/960  m.  The 
horizontal  solid  lines  are  the  steady  state  amplitudes  Ai  i™  and  A2)oTH  from  the 
theory  of  Chaturvedi  and  Ossakow  [1979]. 
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OICY  ATTN  D.  CLARK 

MARTIN  MARIETTA  CORP 
ORLANDO  DIVISION 
P.  O.  BOX  5837 
ORLANDO,  FL  32805 

OICY  ATTN  R.  HEFFNER 

MCDONNELL  DOUGLAS  CORPORATION 
5301  BOLSA  AVENUE 
HUNTINGTON  BEACH,  CA  92647 
OICY  ATTN  N.  HARRIS 
OICY  ATTN  J.  MOULS 
OICY  ATTN  GEORGE  MROZ 
OICY  ATTN  W.  OLSCN 
OICY  ATTN  R.  W.  MALPRIN 
OICY  ATTN  TECHNICAL  LIBRARY  SERVICES 


MISSION  RESEARCH  CORPORATION 
735  STATE  STREET 


SANTA  BARBARA,  CA  93101 

OICY 

ATTN  P.  FISCHER 

OICY 

ATTN  W.  F.  CREVIER 

OICY 

ATTN  STEVEN  L.  GUTSCHE 

OICY 

ATTN  D.  SAPPENFIELO 

OICY 

ATTN  R.  BOGUSCH 

OICY 

ATTN  R.  HENDRICK 

OICY 

ATTN  RALPH  KILB 

OICY 

ATTN  DAVE  SOWLE 

OICY 

ATTN  F.  FAJEN 

OICY 

ATTN  M.  SCHEI6E 

OICY 

ATTN  CONRAD  L.  LONGMIRE 

OICY 

ATTN  WARREN  A.  SCHLUETER 

MITRE  CORPORATION,  THE 

P.  0.  BOX 

208 

BEDFORD,  MA  01730 

OICY 

ATTN  JOHN  MORGANSTERN 

OICY 

ATTN  G.  HAROING 

OICY 

ATTN  C.  E.  CALLAHAN 

JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MD  20810 

OICY  ATTN  DOCUMENT  LIBRARIAN 
OICY  ATTN  THOMAS  POTEMRA 
OICY  ATTN  JOHN  DASSOULAS 

LOCKHEEO  MISSILES  &  SPACE  CO  INC 
P.  O.  BOX  504 
SUNNYVALE,  CA  94088 

OICY  ATTN  OEPT  60-12 
OICY  ATTN  D.  R.  CHURCHILL 

LOCKHEEO  MISSILES  AND  SPACE  CO  INC 
3251  HANOVER  STREET 
PALO  ALTO,  CA  94304 

OICY  ATTN  MARTIN  WALT  DEPT  52-10 
OICY  ATTN  RICHARO  G.  JOHNSON  OEPT  52-12 
OICY  ATTN  W.  L.  IWO F  OEPT  52-12 


MITRE  CORP 

WESTGATE  RESEARCH  PARK 
1820  OOLLY  MAO I  SON  BLVD 
MCLEAN,  VA  22101 

OICY  ATTN  W.  HALL 
OICY  ATTN  W.  FOSTER 

PACIFIC-SIERRA  RESEARCH  CORP 
1456  CLOVERF IELO  BLVD. 

SANTA  MONICA,  CA  90404 

OICY  ATTN  E.  C.  F IELO  JR 

PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 

(NO  CLASS  IF IEO  TO  THIS  ADDRESS) 

OICY  ATTN  IONOSPHERIC  RESEARCH  LAB 
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COMAKER 

NAVAL  OCEAN  SYSTEMS  CENTER 
SAN  01  EGO,  CA  92152 

03CY  ATTN  CODE  532  W.  MOLER 
01CY  ATTN  CODE  0230  C.  BAGGETT 
01CY  ATTN  COOE  81  R.  EASTMAN 

DIRECTOR 

NAVAL  RESEARCH  LABORATORY 
WASHINGTON,  D.C .  20375 

01CY  ATTN  CODE  6700  TIMOTHY  P.  COFFEY 

(25  CYS  IF  UNCLASS,  1  CY  IF  CLASS) 

0 ICY  ATTN  CODE  6701  JACK  0.  BROWN 
0 ICY  ATTN  CODE  6780  BRANCH  HEAD  (150  CYS 

IF  UNCLASS,  1  CY  IF  CLASS) 

0 ICY  ATTN  CODE  7500  HQ  COMM  DIR  BRUCE  WALO 

0 ICY  ATTN  CODE  7550  J.  OAV1S 

0 ICY  ATTN  CODE  7580 

01CY  ATTN  COOE  7551 

0 ICY  ATTN  CODE  7555 

0 ICY  ATTN  CODE  67  30  E.  MCLEAN 

01CY  ATTN  COOE  7127  C.  JOHNSON 

COMMANDER 

NAVAL  SEA  systems  COMMAND 
WASHINGTON,  O.C.  20362 

0 ICY  ATTN  CAPT  R.  PITKIN 

COMMANDER 

NAVAL  SPACE  SURVEILLANCE  SYSTEM 
DAHLGREN,  VA  22448 

01CY  ATTN  CAPT  J.  H,  BURTON 

OFFICER- IN-CHARGE 
NAVAL  SURFACE  WEAPONS  CENTER 
WHITE  OAK,  SILVER  SPRING,  MD  20910 
01CY  ATTN  CODE  FJ1 

01  RECTOR 

STRATEGIC  SYSTEMS  PROJECT  OFFICE 
oepartment  OF  THE  NAVY 
WASHINGTON,  D.C.  20376 
0 ICY  ATTN  NSP-2141 
01CY  ATTN  NSSP-2722  FRED  WIMBERLY 

NAVAL  SPACE  SYSTEM  ACTIVITY 
P.  0.  BOX  92960 
WORLDWAY  POSTAL  CENTER 
LOS  ANGELES,  CALIF.  90009 

0 ICY  ATTN  A.  B.  HAZZARO 

HEADQUARTERS 

ELECTRONIC  SYSTEMS  DIVISION/OC 
DEPARTMENT  OF  THE  AIR  FORCE 
HANSCOM  AFB,  MA  01731 

0 ICY  ATTN  OCKC  MAJ  J.  C.  CLARK 


COMMANDER 

FOREIGN  TECHNOLOGY  DIVISION,  AFSC 
WRIGHT-PATTERSON  AFB,  OH  45433 
01CY  ATTN  NICO  LIBRARY 

0 ICY  ATTN  ETOP  3.  BALLARO 

COMMANDER 

ROME  AIR  OEVELOPMENT  CENTER,  AFSC 
GRIFF ISS  AFB,  NY  13441 

01CY  ATTN  DOC  L 1BRARY/TSLD 
01CY  ATTN  OCSE  V.  COYNE 

SAMSO/SZ 

POST  OFFICE  BOX  92960 
WORLDLY  POSTAL  CENTER 
LOS  ANGELES,  CA  90009 
(SPACE  DEFENSE  SYSTEMS) 

0 ICY  ATTN  SZJ 


STRATEGIC  AIR  COMMAND/XPFS 
OFFUTT  AFB,  NB  68115 

01CY  ATTN  XPFS  MAJ  B.  STEPHAN 
01CY  ATTN  ADWATE  MAJ  BRUCE  BAUER 
01CY  ATTN  NRT 

01CY  ATTN  DOK  CHIEF  SCIENTIST 

SAMSO/YA 
P.  0.  BOX  92960 
WORLDWAY  POSTAL  CENTER 
LOS  ANGELES,  CA  90009 

01CY  ATTN  YAT  CAPT  L.  BLACKWELDER 

SAMSO/SK 
P.  0.  BOX  92960 
WORLDWAY  POSTAL  CENTER 
LOS  ANGELES,  CA  90009 

01CY  ATTN  SKA  (SPACE  COMM  SYSTEMS)  M.  CLAV I N 
SAMSO/MN 

NORTON  AFB,  CA  92409 
(MINUTEMAN) 

01CY  ATTN  MNNL  LTC  KENNEOY 
COMMANDER 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
HANSCOM  AFB,  MA  01731 

01CY  ATTN  EEP  A.  LORENTZEN 


OEPARTMENT  OF  ENERGY 


OEPARTMENT  OF  ENERGY 
ALBUQUERQUE  OPERATIONS  OFFICE 
P.  0.  BOX  5400 
AL3UQUERQUE,  NM  87115 

01CY  ATTN  DOC  CON  FOR  D.  SHERWOOD 

DEPARTMENT  OF  ENERGY 
LIBRARY  ROOM  G-042 
WASHINGTON,  D.C.  20545 

01CY  ATTN  OOC  CON  FOR  A.  LABOWITZ 

EG&G,  INC. 

LOS  ALAMOS  DIVISION 

P.  0.  BOX  809 

LOS  ALAMOS,  NM  85544 

01CY  ATTN  OOC  CON  FOR  J.  BREEDLOVE 

UNIVERSITY  OF  CALIFORNIA 
LAWRENCE  LIVERMORE  LABORATORY 
P.  0.  BOX  808 
LIVERMORE,  CA  94550 

01CY  ATTN  DOC  CON  FOR  TECH  INFO  OEPT 
01CY  ATTN  DOC  CON  FOR  L-389  R.  OTT 
01CY  ATTN  DOC  CON  FOR  L-31  R.  HAGER 
01CY  ATTN  DOC  CON  FOR  L-46  F.  SEWARD 

LOS  ALAMOS  SCIENTIFIC  LABORATORY 

P.  0.  BOX  1663 

LOS  ALAMOS,  NM  87545 


01CY 

ATTN 

DOC 

CON 

FOR 

J. 

WOLCOTT 

01CY 

ATTN 

DOC 

CON 

FOR 

R. 

F.  TASCHEK 

01CY 

ATTN  DOC 

CON 

FOR 

E. 

JONES 

01CY 

ATTN  DOC 

CON 

FOR 

J. 

MALIK 

01CY 

ATTN 

DOC 

CON 

FOR 

R. 

JEFFRIES 

0  ICY 

ATTN  DOC 

CON 

FOR 

J. 

ZINN 

01CY 

ATTN 

OOC 

CON 

FOR 

P. 

KEATON 

01CY 

AT  TN 

DOC 

CON 

FOR 

D. 

westervelt 

SAND I A  LABORATORIES 

P.  0.  BOX 

5800 

ALBUQUERQUE,  NM 

87115 

01CY 

ATTN 

OOC 

CON 

FOR 

J. 

-MARTIN 

01CY 

ATTN 

OOC 

CON 

FOR 

w. 

BROWN 

01CY 

ATTN 

OOC 

CON 

FOR 

A. 

THORNS  R'JUGh 

01CY 

ATTN 

DOC 

CON 

FOR 

T. 

WRIGHT 

01CY 

ATTN 

OOC 

CON 

FOR 

D. 

oamlgren 

01CT 

ATTN 

OOC 

CON 

FOR 

3141 

01CY 

ATTN 

DOC 

CON 

FOR 

SPACE  pROuECT  DIV 
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SAND I A  LABORATORIES 
LIVERMORE  LABORATORY 
P.  0.  BOX  969 
LIVERMORE,  CA  94550 

01CY  ATTN  l)OC  CON  FOR  3.  MU RP HEY 
01CY  ATTN  OOC  CON  FOR  T.  COOK 

OFFICE  OF  MILITARY  APPLICATION 
OE  PART  ME  NT  OF  ENERGY 
WASHINGTON,  O.C.  20545 

01CY  ATTN  OOC  CON  FOR  0.  GALE 


OTHER  GOVERNMEiNT 


CENTRAL  INTELLIGENCE  AGENCY 
ATTN  RO/SI,  RM  5G48,  HQ  BLDG 
WASHINGTON,  O.C.  20505 

01CY  ATTN  OSI/PSIO  RM  5F  19 

DEPARTMENT  OF  COMMERCE 

national  bureau  of  stakoards 

WASHINGTON,  O.C.  20234 

(ALL  CORRES:  ATTN  SEC  OFFICER  FOR) 
01CY  ATTN  R.  MOORE 


INSTITUTE  FOR  TELECOM  SCIENCES 
NATIONAL  TELECOMMUNICATIONS  $  INFO  ADMIN 
BOULDER,  CO  80303 

01CY  ATTN  A.  JEAN  (UNCLASS  ONLY) 
01CY  ATTN  W.  UTLAUT 
0 ICY  ATTN  D.  CROMBIE 
01CY  ATTN  L.  BERRY 

NATIONAL  OCEANIC  &  ATMOSPHERIC  ADMIN 
ENVIRONMENTAL  RESEARCH  LABORATORIES 
DEPARTMENT  Of  COMMERCE 
BOULDER,  CO  80302 

0 ICY  ATTN  R.  GRUBB 

Q1CY  ATTN  AERONOMY  LAB  G.  REID 


AEROSPACE  CORPORATION 


).  BOX  92957 
ANGELES,  CA 

90009 

OlCY 

ATTN 

I .  GARFUNKEL 

OlCY 

ATTN 

T.  SALMI 

OlCY 

ATTN 

V.  JOSEPHSON 

OlCY 

ATTN 

S.  BOWER 

OlCY 

ATTN  N.  STOCKWELL 

OlCY 

ATTN  0.  OLSEN 

OlCY 

ATTN 

J.  carter 

OlCY 

ATTN 

F.  MORSE 

OlCY 

ATTN 

SMFA  FOR  PWW 

ANALYTICAL  systems  ENGINEERING  CORP 
5  OLO  CONCORD  ROAD 
BURLINGTON,  MA  01803 

01CY  ATTN  RAO IO  SCIENCES 

BERKELEY  RESEARCH  ASSOCIATES,  INC. 
P.  O.  BOX  983 
BERKELEY,  CA  94701 

01CY  ATTN  J.  'WORKMAN 


BOEING  COMPANY,  THE 
P.  0.  BOX  3707 
SEATTLE,  wA  98124 

01CY  ATTN  G.  KEISTER 
01 CY  ATTN  D.  MURRAY 
01CY  ATTN  G.  HALL 
01CY  ATTN  J.  KENNEY 

CALIFORNIA  AT  SAN  DIEGO,  UN IV  OF 
IPAPS,  B— 0 1 9 
LA  JOLLA,  CA  92093 

01CY  ATTN  HENRY  G.  BOOKER 

BROWN  ENGINEERING  COMPANY,  INC. 

CUMMINGS  RESEARCH  PARK 
HUNTSVILLE,  AL  35807 

01CY  ATTN  ROMEO  A.  OELIBERIS 

CHARLES  STARK  DRAPER  LABORATORY,  INC. 

555  TECHNOLOGY  SQUARE 
CAMBRIDGE,  MA  02139 

OlCY  ATTN  D.  8.  COX 
OlCY  ATTN  J.  P.  GILMORE 

COMPUTER  SCIENCES  CORPORATION 
5565  ARLINGTON  3LVD 
FALLS  CHURCH,  VA  22046 
OlCY  ATTN  H.  BLANK 
OlCY  ATTN  JOHN  SPOOR 
OlCY  ATTN  C.  NAIL 

COMSAT  LABORATORIES 
LI NTH  I  CUM  ROAD 
CLARKSBURG,  MD  20734 

OlCY  ATTN  G.  HYDE 

CORNELL  UNIVERSITY 

DEPARTMENT  OF  ELECTRICAL  ENGINEERING 
ITHACA,  NY  14850 

OlCY  ATTN  D.  T.  FARLEY  JR 

ELECTROSPACE  SYSTEMS,  INC. 

BOX  1359 

RICHARDSON,  TX  75Q&0 

OlCY  ATTN  H.  LOGSTON 

OlCY  ATTN  SECURITY  (PAUL  PHILLIPS) 

ESl  INC. 

495  JAVA  DRIVE 
Sunnyvale,  ca  94086 

OlCY  ATTN  J.  ROBERTS 
OlCY  ATTN  JAMES  MARSHALL 
OlCY  ATTN  C.  W.  PRETTIE 

FORD  AEROSPACE  S  COMMON I  CAT  I ONS  CORP 
3939  FABIAN  WAY 
PALO  ALTO,  CA  94303 

OlCY  ATTN  J.  T.  MATTINGLEY 

GENERAL  ELECTRIC  COMPANY 
SPACE  DIVISION 
VALLEY  FORGE  SPACE  CENTER 
GODOARD  BLVD  KING  OF  PRUSSIA 
P.  0.  BOX  8555 
PHILADELPHIA,  PA  19101 

OlCY  ATTN  M.  H.  BORTNER  SPACE  SCI  LAB 

GENERAL  ELECTRIC  COMPANY 
P.  0.  BOX  1122 
SYRACUSE,  NY  13201 

OlCY  ATTN  F.  RE I3ERT 

GENERAL  ELECTRIC  COMPANY 
TEMPO-CENTER  FOR  ADVANCED  STUDIES 
316  STATE  STREET  (P.O.  DRAWER  QQ) 

SANTA  BARBARA,  CA  93102 
OlCY  ATTN  DAS  I  AC 
OlCY  ATTN  DON  CHANDLER 
OlCY  ATTN  TOM  BARRETT 
OlCY  ATTN  TIM  STEPHANS 
OlCY  ATTN  WARREN  S.  KNAPP 
OlCY  ATTN  WILLIAM  MCNAMARA 
OlCY  ATTN  3.  GAMS  ILL 
OlCY  ATTN  MACK  STANTON 
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PhOTOMcTRICS,  INC. 

44  2  MARR-TT  ROAD 
LEXINGTON,  MA  32173 

01CY  ATTN  IRVING  <Cr SKY 

PHYSICAL  DY\AMI  CS  INC. 

P.  G.  30 X  3027 
BE-LtVUt,  WA  93009 

0 ICY  ATTN  c.  J.  FrSMOjW 


iCIENCc.  AP-L  I A7  |.',NS ,  I  VlDRPCRA ' 
3400  jRI,t 

mC-EAN,  /A  2 2101 

01  _Y  AT "N  j.  C DGK A f \£ 

SCIENCE  APPLICATIONS,  INC. 

30  .MISSION  JRIv  = 

PLEASANTON,  CA  s4S>ou 
01CY  ATTN  SZ 


PHYSICAL  DYNAMICS  INC. 
p.  O.  SOX  1069 
BERKELEY,  CA  94701 

01CY  ATTN  A.  T-QMPSCN 

PSD  ASSOCIATES 

P.  0.  box  9695 

MARINA  DEL  REY,  CA  90291 

0 ICY  ATTN  FORREST  oILMORE 
01CY  ATTN  BRYAN  GABBARD 
01 CY  ATTN  WILLIAM  5.  WRIGHT  JR 
01CY  ATTN  ROBERT  F.  lELEVIER 
01CY  ATTN  WILLIAM  J.  < AREAS 
0 ICY  ATTN  H.  ORY 
01CY  ATTN  C.  MACDONALD 
01CY  ATTN  R.  TURCO 

rand  corporation,  THE 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 

0 ICY  ATTN  CULLEN  CRAIN 
01CY  ATTN  ED  BEDROZ I  AN 

RIVERSIDE  RESEARCH  INSTITUTE 
80  WEST  END  AVENUE 
NEW  YORK,  NY  10023 

0 ICY  ATTN  VINCE  TRAPANI 

SCIENCE  APPLICATIONS,  INC. 

P.  O.  BOX  2351 
LA  JOLLA,  CA  92038 

01 CY  ATTN  LEWIS  M.  LINSCN 
0 ICY  ATTN  DANIEL  A.  HAMLIN 
01CY  ATTN  0.  SACHS 
0 ICY  ATTN  E.  A.  STRAKER 
01CY  ATTN  CURTIS  A.  SMITH 
01 CY  ATTN  JACK  MCDOUGALL 

RAYTHEON  CO. 

528  BOSTON  POST  ROAO 
SUDBURY,  MA  01776 

01CY  ATTN  BARBARA  ADAMS 


SRI  INTERNATIONAL 
333  RAVENS ^000  AVENUE 
MENLO  PARK,  CA  9*325 

01CY  ATTN  DONA^j  NEILSON 
Q1CY  ATTN  ALAN  BURNS 
0 ICY  ATTN  G.  SMITH 
01CY  ATTN  L.  L.  COBB 
dlCY  ATTN  DAVID  A.  JOHNSON 
0 ICY  ATTN  .vALTER  G.  ChtSNUT 
01CY  ATTN  CHARLES  L.  RING 
01CY  ATTN  WALTER  JAYE 
01CY  ATTN  M.  BARON 
01CY  ATTN  RAY  L.  LEADABRAND 
0 ICY  ATTN  G.  CARPENTER 
0 ICY  ATTN  G.  PRICE 
0 ICY  ATTN  J.  PETERSON 
01CY  ATTN  R.  HAKE,  JR. 

01CY  ATTN  V.  GONZALES 
01CY  ATTN  D.  MCDANIEL 

TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
BEDFORD,  MA  01730 

01CY  ATTN  W.  P.  BOQUIST 

TRW  DEFENSE  6  SPACE  SYS  GRDUP 

ONE  SPACE  l=ARK 

REDONDO  BEACH,  CA  90278 

01CY  ATTN  R.  K.  PLEBUCH 
01CY  ATTN  S.  ALTSCHULER 
01CY  ATTN  D.  OEE 

VISIDYNE,  INC. 

19  THIRD  AVENUE 

NORTH  WEST  INDUSTRIAL  PARK 

BURLINGTON,  MA  01803 

01CY  ATTN  CHARLES  HUMPHREY 
01CY  ATTN  J.  W.  CARPENTER 


SCIENCE  APPLICATIONS,  INC. 
HUNTSVILLE  DIVISION 
2109  W.  CLINTON  AVENUE 
SUITE  700 

HUNTSVILLE,  AL  35805 

0 ICY  ATTN  DALE  H.  DIVIS 
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